Exact calculations of first-passage quantities on recursive networks 
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We present general methods to exactly calculate mean-first passage quantities on self-similar 
networks defined recursively. In particular, we calculate the mean first-passage time and the splitting 
probabilities associated to a source and one or several targets; averaged quantities over a given set 
of sources (e.g., same-connectivity nodes) are also derived. The exact estimate of such quantities 
highlights the dependency of first-passage processes with respect to the source-target distance, which 
has recently revealed to be a key parameter to characterize transport in complex media. We explicitly 
perform calculations for different classes of recursive networks (finitely ramified fractals, scale-free 
(trans) fractals, non- fractals, mixtures between fractals and non-fractals, non-decimable hierarchical 
graphs) of arbitrary size. Our approach unifies and significantly extends the available results in the 
field. 



As they account for transport efficiency in complex 
media, first-passage processes [T] on fractal or complex 
networks have given rise to a lot of interest in the past 
few years. Such networks, showing a very broad spec- 
trum of topological and thus transport characteristics, 
have been opportunely used to model a wide range of 
systems, from biology to sociology or computer science 
[2]-[5]. Within this outline, the exact determination of 
first-passage quantities has been the topic of many works 
[6HT3], in a dialog with parallel investigations of general 
scalings and asymptotic results [141117^ . The global mean 
first-passage time (GMFPT), defined as the mean first- 
passage time of a random walk starting from an arbitrary 
site (source) in presence of a fixed trap (target), and in 
particular its scaling with the number of nodes was 
meant to encompass the general properties of transport 
on fractal and transfractal networks. Yet, recent results 
[16 have highlighted the insufficiency of this mere de- 
scription and put forward the role of non-averaged quan- 
tities, associated to a single starting point. For example, 
in the context of diffusion limited reactions in complex 
media, the initial position of the reactants has indeed 
been shown to be a key parameter that can control the 
entire kinetics of the process [18 . 

In that context, we present a general algebraic method 
to calculate exactly first-passage quantities on self-similar 
networks, for a given source point. More precisely, we 
consider here two first-passage observables: 

• the mean first-passage time (MFPT) from a site S 
to a site T, or to several sites {T^}, i.e. the average 
time it takes a random walker that starts at S to 
reach T or any point in {T^}; 

• the splitting probability, defined in presence of sev- 
eral targets Ti as the probability, starting from 
to reach a given target T^^ before all the other tar- 
gets. 

Actually, these quantities both satisfy simple but formal 
N X N linear systems, where N is the number of sites of 



the network. The regime of interest is typically N large 
and makes the explicit resolution of such linear systems 
out of reach. The aim of the present work is to provide a 
general method, applicable to a broad class of networks, 
which yields explicit and exact formulas for the MFPT 
and splitting probabilities even for very large N. To do 
so, we make use of the self-similar properties of the net- 
works considered and develop a general renormalization 
scheme following ideas presented in ^llj . 

More precisely, we consider in this paper self-similar 
networks, which can be defined recursively ^ by consid- 
ering the result of the action of a similarity transforma- 
tion r(r) that maps all points i of a network F onto new 
points = r{i). Considering that the transformation r is 
homothetic of ratio p < 1, a network is called self-similar 
if at generation n its iteration F^ is equal to the union of 
p replicas of r(F^_i). As a result, self-similar networks 
can be split into a finite number of equal sub-units, those 
sub-units being the network at the previous generation. 
Explicit examples will be given throughout this paper. 

An important point is the relationship between self- 
similarity and fractality. A network is called fractal if 
one can define a constant df such as A/" ~ R^f ^ where 
N denotes the volume (or number of nodes) and R the 
chemical diameter. Note that whereas most examples 
of deterministic fractal networks are self-similar, the re- 
ciprocal is false: there exists self-similar networks that 
do not show the fractal property, such as the so-called 
(-u, 'u) -flower networks with u = 1 that will be defined 
more precisely below. In the latter example, the diameter 
scales as a logarithm of the number of nodes: R ~ log N ^ 
which is referred to as the small-world property. Such a 
network can be formally seen as a fractal network with 
infinite fractal dimension, and is sometimes called trans- 
finite fractal or transfractal [19 . 

Beyond fractality, self-similar networks may also ex- 
hibit other prominent properties that seem to be com- 
mon to real networks, especially biological and social 
networks, such as scale-free features or modular struc- 
ture [20^. The former implies that the distribution of 
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FIG. 1: Examples of self similar networks : 4tli generation 
Sierpinski gasket and 5tli generation T-graph. 



the degree of the nodes of the network follows a power 
law, while the latter means that the network can be 
divided into groups (modules), within which nodes are 
more tightly connected with each other than with outside 
nodes [21^.22^. In the last part of this article we focus on 
a class of hierarchical, non-decimable, recursive networks 
able to capture simultaneously scale-free behavior and 
modular structure, yet preserving (weak) self-similarity 
[201 [21]. As we will see, the calculation of first-passage 
quantities on such networks requires an alternative ap- 
proach that will be presented in the last section of this 
article. 

This paper is organized as follows. In the first Section, 
we introduce general definitions and present in detail 
the method of calculation of splitting probabilities and 
MFPTs on the example of the Sierpinski gasket. Then, 
we extend this method to other deterministic self-similar 
networks on the example of the Song-Havling-Makse net- 
works (Section II), while further examples (T graph and 
{u^ v) flowers) are given in Appendix. We stress that this 
approach allows one to calculate explicitly the splitting 
probabilities and MFPTs for any starting site of the lat- 
tice (but for specific targets). In addition, we will show 
that it leads to simple expressions of selective averages 
over starting sites (that is averages over specific class of 
starting sites, to be defined below). In the particular case 
where the average is performed over all starting sites, we 
recover expressions of global MFPTs recently obtained. 
Finally, in Section III, we present an alternative method 
of calculation in the case of hierarchical networks. 



I. THE GENERAL METHOD ILLUSTRATED 
ON THE EXAMPLE OF THE SIERPINSKI 
GASKET 

A. Definitions 

In this section, we first introduce the class of networks 
to which our method of calculation can be applied, and 
give basic definitions that will be used throughout the 
paper. Note that here we do not aim at giving mathe- 
matically formal definitions, and will rather largely rely 
on explicit examples. 

Hierarchical networks. As stated in introduction, we 
consider hierarchical networks which can be defined [2] by 
considering the action of a rescaling transformation r(r) 
that maps all points z of a network F onto new points 
i^ = r{i). A self-similar network is then constructed 
recursively from an elementary motif (initiator) Fq by 
writing F^ = Ui=i..pr(F^_i). For example, the Sierpin- 
ski gasket of generation n is built by joining p = 3 copies, 
called subunits, of Sierpinski gaskets of generation n — 1 
(see Figs. [l]and[2|, where the initiator is the elementary 
triangular network with three nodes. 




FIG. 2: Sierpinski gasket : renormalization scheme and asso- 
ciated crossing times Tk. 

Levels. A given site io of a hierarchical network F^ is 
said to belong to the level ^, if there exists an integer ^ < 
n such that io ^ P^. We will denote by Ck the set of nodes 
of level k. Note that the networks are constructed such 
that a site that belongs to the level ^ also belongs to levels 
^ + 1, ^ + 2, n . Figure [3] illustrates this definition 
for the 3rd-generation of a Sierpinski gasket (see legend) . 
Note that, for the nth-generation of a Sierpinski gasket, 
there are - — sites on the level k {k G 1, ...,n). 

Labeling of the sites and subunits. Any subunit of 
level /c in a self-similar network of generation n can be 
reached recursively by defining a path^ i.e. a sequence 
{zo,...,^/c} where ij labels the position of each of the 
nodes of the initiator Fq. An example is given in Fig- 
ure [4] in the case of the Sierpinski gasket, by assigning 
the value i^ = to a top subunit, i/e = 1 to a left one 
and i/c = 2 to a right one. With these rules, the path 
{^17^27^37^} = {2,0,1,2} allows to locate the last-level 
sub-unit {t^^\l'^^\r'^^^) within the 4th level. 

Transport process. We consider throughout this pa- 
per a nearest neighbor Markovian random walker char- 
acterized by generic transition probabilities w{r'\r) from 
site r to r^ Unless specified we will consider isotropic 
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FIG. 3: Color online. Third-generation Sierpinski gasket : 
levels. Blue circles : level 3 ; green squares: levels 2,3 ; ma- 
genta stars: levels 1,2,3 ; orange diamonds: levels 0,1,2,3. 



^(0) ^(2) 




ii = 2(=right) 
\2 = 0(= top) 

13 = l(=left) 

14 = 2 (= right) 



FIG. 4: Color online. Sierpinski gasket : example of labels. 



random walks such that w{y'\y) = l//^(r) where /^(r) is 
the connectivity of node r. 

B. Splitting probabilities 

In this section we wish to calculate the splitting prob- 
ability Pyt\va (^)' defined as the probability for a Marko- 
vian random walker starting at r to reach the target 
in the presence of other absorbing sites {r^^.j^, in other 
words the probability to reach the site before all the 
other absorbing sites. For the sake of readability, we will 
make use of the following notation: 

^'rHr.,(r)=P(r). 

We will present the method on the example of the Sier- 
pinski gasket and show later on how it can be generalized 
to any self-similar network. Let us consider a general sub- 
unit A/e_i at a given level /c — 1, which is depicted in the 
right hand side of Fig. [4j We assume that the splitting 
probability P(r) is known for all starting sites of level 
k — 1 in A/e_i (in the Sierpinski gasket there are only 3 
such sites which correspond to the summits A^B^C of 
the main triangle defining A/c_i), and that the absorbing 



sites yt and {r^^ji are located outside A/e_i. Here the 
subunit A/c_i is the union of p = 3 copies of subunits A^. 
The splitting probabilities for sites r of level k in A/^-i 
(namely the sites A, B, C, D, E, F in Fig |2| satisfy the 
following backward equation pLj: 

0= Hr'\r)P(r')-P{r). (1) 

Here 7r(r'|r) is the splitting probability that the random 
walker starting from the node r of level k in A^-i reaches 
first the site among all sites of level k in A/c_i. The 
sites of level k in Ak-i actually form a graph Fi of gen- 
eration 1 up to a rescaling factor (see Fig|2|. In the case 
of an isotropic walk, it is clearly seen on the example 
of the Sierpinski gasket (for which all nodes of the ini- 
tiator Fo are equivalent) that 7r(r'|r) = K;(r'|r), where 
w{r'\r) is the elementary transition probability on Fi. 
More explicitly, on Fig |2) on has 7t{B\E) = 7t{D\E) = 
7r{F\E) = 7r{A\E) = 1/4. In the remainder of the arti- 
cle, we consider isotropic random walks on networks, for 
which this property holds. Nevertheless, the calculation 
method that we present can be in principle extended to 
directed or non-uniformly weighted networks (such that 
7r(r'|r) ^ 7r(r'|r)), as long as the scale-invariance hypoth- 
esis is fulfilled. 

In this example, let us assume that P(A), P{B) and 
P{C) are known; the expressions of P{D), P{P^) 
P{F) of the splitting probabilities starting from the 
points D^E^F of level k in A^-i can then be obtained 
readily by making use of Eq. ([T]). One has 

P{E)=^-[P{A) + P{F) + P{D) + P{B)], (2) 
P{F)=^-[P{A) + P{E) + P{D) + P{C)], (3) 
P{D) = i [P{B) + P{E) + P{F) + P{C)] , (4) 

which can be rewritten as a linear system 

P{D)\ /1/5 2/5 2/5\ (P{A)\ 
P{E) = 2/5 2/5 1/5 P{B) . (5) 
P{F)J \2/5 1/5 2/5/ \P{C)J 

We now proceed iteratively and consider for example 
the upper subunit of Aj^-i that contains the nodes 
A, F of level k. This choice is taken into account by 
assigning to a variable ik the value = top. It is of course 
possible to apply the same operation to any of the two 
other subunits : the lower-left one {ik = I = left) or the 
lower-right one {i^ = 2 = right). Let us then rename the 
nodes of level /c — 1 in Ak-i and the nodes of level k in 
A/c as follows: 

^(/c-l) ^ iik-1) ^ ^(/c-1) ^ 
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Equation ([5| can then be rewritten as: 




1 

2/5 2/5 1/5 
2/5 1/5 2/5 




(/c-l) 



(6) 



Similarly, one can define three matrices M^i^ such that 
equation (|6| reads for each value of ik corresponding to 
either a top, left or right subunit: 



(7) 



with 
Mo 



1 

2/5 2/5 1/5 
,2/5 1/5 2/5. 



and M2 



Ml 



^2/5 2/5 1/5^ 
1 
1/5 2/5 2/5. 



^2/5 1/5 2/5^ 
1/5 2/5 2/5 
1 



(8) 



As shown above, any subunit (of level k) in the net- 
work can be reached recursively by defining a path^ i.e. 
a sequence with ii G {0,1,2} for 1 < / < /c. 

Iterating equation Q then yields straightforwardly: 



. . . A/f- -pC^o) 



(9) 



This shows that as soon as V^^^^ is known for a given level 
/co, the splitting probability starting from any point of 
level k > ko can be obtained exactly and only requires to 
compute a product of k — ko 3x3 matrices. In particular 
if the targets are chosen among the 3 sites of level then 
7^*^^°^ is calculated trivially and the splitting probability 
for any starting point of the network is readily obtained. 

Let us consider an explicit example illustrated in Fig. 
|4]that represents a Sierpinski gasket of generation 4. We 
aim at calculating the splitting probability to reach t*^^^ 
before r^^^ starting from /^^^: 

P,C0)|H0)aW)^PaW). 

P(/^^^) is the 2nd coordinate of the vector V^^\ associ- 
ated to the subunit of level 4 (with vertices in red) . Since 
the three matrices Mq, Mi and M2 are known, all we 
need to determine is the path of that subunit and the 
value of 7^^^^. As before, one has: 



and clearly 



n,^2,^3,u = 2,0,1,2 




Thus, applying equation (|9|), we get V^^"^ 
/234/625\ 

229/625 , so that finally: 
\ 89/250 / 



229 
625 



(10) 



C. MFPT 

In this section we consider the MFPT of a random 
walker to a target site starting from a site r, that we 
denote T(r). Following the main steps of derivation of 
the splitting probability, we first write down a backward 
equation for the MFPT starting from a given site r of 
level /c in a subunit A/^-i of level k — 1 (see Fig. [2|: 



= AT(r)= ^(r|r)T(rO-T(r). (11) 



As discussed before, the splitting probability 7r(r'|r) is 
readily given by the transition probability w{r'\r) on the 
corresponding graph Fi of generation 1 formed by the 
sites of level k in A^-i (sites A, B, C, D, E, F in Fig[2|. 
In addition we introduced the quantity Tk defined as the 
time it takes a random walker to exit a subunit of level 
k. On the example of the right hand side of Fig|2j Tk is 
the mean time to reach either B^D^F^A starting from 
E. Note that by symmetry it can also be defined as 
the mean time to reach either E or F starting from A. 
In this example, all the exit nodes of a given subunit 
play a symmetric role, and by construction, all subunits 
of a given level are the same. Therefore Tk does not 
depend on r in equation (11). This property results here 



from the symmetry of the initiator Fq, in which all nodes 
are equivalent, and from the symmetry of the transition 
probabilities. A stated previously, we will consider in this 
paper only networks having this property. 

More explicitly, we rely on the example of Fig. [2] and 
assume that the MFPT starting from the sites A, 5, C of 
level /c — 1 is known. Using Eq. (11), one can write : 



1 r 



T{E) = -[(T(A)+r,) + (T(F)+r,) 
+(T(I))+rfc) + (T(5)+rfc)l, 



and the the 2 similar relations at nodes F and D. Fol- 
lowing the derivation of splitting probabilities above, the 
MFPT starting from any site of level k can be expressed 
linearly in terms of the MFPT starting from sites of level 
k — 1. For example, focusing on the top subunit AEF 
(corresponding to the choice ik = 0), one obtains: 




TkVo^MoT^^ 



1 

2/5 2/5 1/5 

2/5 1/5 2/5y 

1) 



(12) 



(nt)\ 

where T^^^ = T{1) denotes the vector of MFPTs 
\T{r)J 

starting from the 3 vertices of level A: of a given subunit of 
level k located by its path {zi, ...,i/c}. Similar equations 
can be obtained for the left {ik = 1) or the right {ik = 2) 
subunit, and yield the following general recursive relation 
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T'^"^ =TkV,,+M,,T^''-'\ (13) 
where the Aiif^ matrices are given by Eq. (|8| and 



Vo 



Vi 



(14) 



Assuming that 7"^^°^ is known for some level /cq, the 
MFPT for a given level k > ko can then be written as: 



Example 2 : a class of sources. In the present ex- 
ample, the aim is to study the MFPT dependence with 
respect to the target-source distance. The target is left 
at t^^\ and the class of sources is the sites located on 
the adjacent edge [t^^\l^^^].^ at the distances = 2^, 
< p < n. Those sources are the lower-left vertices of 
the triangles corresponding to the paths {0} (r = 2^), 
{0,0} (r = 2^-1), {0,0,0} (r = 2^-^), etc. Thus, the 
problem is solved by using formula ([l5| ior k = n — p 
and Mil = -^o for ah We find finally: 



k-l 



r« = Tfe Vi, + Y,n {Mi,Mi,_, ■ ■ • Mi,^,) Vi,, (15) 



l=ko 



where here by definition r^^^Vi^^ = 7"^^°^. We give be- 
low the example where the target is located at the apex, 
for which the determination of T^^^ is straightforward. 
As we show in the next paragraph, can be calculated 
explicitly. In the case where the target is located at the 
apex, Eq. (15) therefore provides an explicit and exact 



expression of the MFPT starting from any starting site 
of the network. Examples will be given below. 

Determination of the exit time Tk • The notations we 
refer to are those of figure |2j It is easily seen that Tk and 
Tk-i are related by the following system: 



Tk-l=Tk + \Cl{E)+1{F)) 

%{E) = i [(Tfe + Tfe_i) + (Tfe + 1(F)) + (Tfe + 1(D)) + Tfe] 
1(D) = I [{Tk + 1{E)) + {Tk + 1{E))] + \Tk 



%{E) = 1(F), 



where and ^{D) denote the mean time to 

reach any of the two B or C sites starting from F 
and D. The solution of this system is r^-i = 5r/e which 
leads to Tk = b'^~^Tn' At level n the subunit is a simple 
triangular graph so that obviously = 1. Finally, one 
has: 



Tk = 5^- 



(16) 



Example 1. Here, we refer again to Fig. [4] and we 
wish to calculate T{1^^^) when the target is set at the 
apex t^^\ First, it is necessary to determine T^^^: 

T(Z(0)) =To+\ T{r(''^) and T(|(o)) = T(r(°)), 
therefore 



T{r = 2P) = 5P (3"-f+i - 1). 



(19) 



The latter result can be compared to the general 
asymptotic expression of the MFPT that has recently 
been derived in in the large system size limit. In 
the case of a target located at the apex of the Sierpinsky 
gasket it reads for N large : 



2iVr^ 



(20) 



where r is the source-target distance, df the fractal di- 
mension of the network and the walk dimension, char- 
acterizing the power-law behaviour of the mean-square 
displacement with respect to time : (Ar^) ~ ^^/d^ 
the present case, 27V = 3^+^ + 3, = ln(3)/ln(2) and 
d^ = ln(5)/ln(2). Thus: 



5^ 3 



n— p+l 



(21) 



As expected this asymptotic regime is recovered by tak- 
ing the large volume limit n ^ oo in the exact expression 



D. Averages 

In this section, we aim at calculating the MFPT to 
a target site averaged over different classes of starting 
points. The average can cover either a class of starting 
points (all the sites of a given level, all sites of a given 
connectivity...) or all the sites of network. In this latter 
case the averaged MFPT is often called the Global MFPT 
[HI [11] . We will show in this section that averages of the 
MFPT over all sites of a given level take simple explicit 
forms. 

In the case of the Sierpinski gasket, we start from Eq. 



( 15 ) that gives an explicit expression of the MFPT start- 



The path and the value of Tk are already known. Finally 
one obtains: 



ing from any of the three points of level /c in a given 
subunit of level k. Since each subunit of level k is in one 
to one correspondence with a path {ii, • • • , ^fe}, one has 
(17) to calculate ... i^}^^^^ - Here we assume that the 

target site is at level and that T^^^ is known. Making 
use of the following identity 



T^^^ = 1176 



(18) 



E E = E E 
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we are back to calculate the expression 

The variable ik is chosen out of a set of p values (in the 
case of the Sierpinski gasket, p = 3), and one has: 

{iu-,ik} 

X {Mo + A^i + ... + Mp)^-' (Vo + Vi + ... + Vp) . 

(22) 

+ Mp and Vtot = Vo + 



We then define M 



tot 



Mo 



... + Vp, and obtain the following general formula: 

k 

paths 1=1 

(23) 

We give below explicit examples. 

Case of a single target at the apex. Let us consider a 
ffiven level k, and denote by 6'^^) the MFPT to a target 
site set at the apex summed over all the starting sites 
of level k. Using expressions (|8| and (14), and the fact 
that p = 3 for the Sierpinski gasket, the diagonalization 
of Mtot leads to: 



J-i 



k-l 





1 




k-l 




(; 


3 




1 1 


B 




1 


3^ 







5fc- 



4 X 3^=-^ I 1 



(24) 



Then, using T'-°'> = tqVq, Eq. Q yields: 
7-(fe) = 4 X 6'^-^ 5"-'^ 1 

paths \ 1 



k-l 




■ (25) 



The first coordinate of the latter expression corre- 
sponds to the sum of the MFPTs starting from the top 
vertex of each of the subunits of level k. Its 2nd and 
3rd coordinates are respectively the sum for the left ver- 
tices and for the right vertices. An examination of the 
Sierpinski gasket shows that the sum of the three co- 
ordinates of the vector defined by Eq. (25) is equal to 
2S'(^) - T(/^o)) _ Knowing that the number of 

sites of level k is we finally obtain the MFPT 



(T) 



averaged over all starting sites of level k: 
ik) 2 



(T) 



3/C+1 ^ I 

5n+i 3^ ^ 4 X 5^ 



-k ok 



3fe+i _^ 1 



(26) 



For the particular case k = n this quantity is the GMFPT 
and is in agreement with the results of [8j. 

The 3 vertices of level are targets. We now assume 
that the 3 vertices of level are targets and wish to calcu- 
late the MFPT averaged over all starting points of level 
k. The only quantity that we need to modify in the cal- 
culation of the previous paragraph is T^^^ . In the case of 
3 targets one has straightforwardly 



7^(0) 



(27) 



The other terms of equation (23), that are given by equa- 



tion (24), are unchanged. We thus have: 



7-(fc) = 5^ 3^-1 _ 5-/c>j ^ I ^ (28) 

paths \ 1 , 



We finally obtain the desired quantity: 
-ik) 5^3^-^1-5"^) 



{Ty 



3/C+1 



E. Conclusion 



(29) 



In this section we have derived on the example of the 
Sierpinski gasket exact expressions of the splitting prob- 
abilities and MFPTs. Our method yields simple expres- 
sions in the case where the target (s) belong(s) to the level 
of the network and applies to any starting site. This 
method can be readily generalized to other self-similar 
networks for which a similar addressing of the subunits 
of each level can be defined. The only task is then to cal- 
culate the matrices A4i^ the vectors Vi and the exit time 
Tfc. The example of the T-graph is given in appendix, 
and further examples that require slight adaptations of 
the method are detailed below. 



II. THE SONG-HAVLIN-MAKSE NETWORKS 

We now consider the example of the Song-Havlin- 
Makse networks, which have been introduced in [23l [24] 
in order to build self-similar networks that mix fractal 
and non-fractal growing schemes. Here, we focus on the 
deterministic version of these networks, and show that 
the method developed in the previous section applies 
upon minor modifications. 

The building scheme of these networks starts from a 
single link between two nodes. The next generation is 
obtained by attaching m new nodes (sons) to each of 
them, then by deleting the original link and by creat- 
ing X new links connecting the newly created sons, as 
shown on figure [5j These networks are self-similar and 
a labeling of each subunit of level k can be defined. As 
before we will derive the corresponding matrices A4if^, 
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FIG. 5: Song-Havlin-Makse network renormalization scheme: 
starting from a link between two nodes, the next generation 
is obtained by attaching m new nodes (sons) to each of them, 
then by deleting the original link and by creating x new links 
connecting the newly created sons. Example of m = 3 and 
X = 1 or X = 2. 



a(k-i) ^{k-D 



i 



FIG. 6: Song-Havlin-Makse network : indexes. The indexa- 
tion of the points does not depend on m nor x (when x / 1), 
because all the x connection links play a symmetric role, like- 
wise the m — X free links. Note the switch of positions of a 
and b in the case ik — 2. 



the vectors V^^ and the exit time r^. It is then possible 
to calculate explicitly splitting probabilities and MFPTs 
for any starting point on the network; here we will calcu- 
late examples of MFPTs averaged over different classes 
of sources. 

Labeling of subunits and sites. Figure |6] shows a possi- 
ble labeling based on 2-dimensional vectors. The scheme 
does not depend on m and x (provided that x > 1). This 
is due to the symmetric role of the x connecting links (re- 
spectively of the m — X free links). Although the results 
do depend on m and this dependence does not apply 
to the matrices and vectors. This can be understood as 
follows. By definition (see Eq. ([7|)), Al^ determines how 
the splitting probabilities P(a*^^^ and P(6'^^^) at level k 
are related to the same quantities at level k — 1. It is easy 
to see that all trajectories starting from one of the nodes 
at level k in one branch will reach sites a^^~^^ or 
of level k — 1 before any other site of other branches. As 
a consequence, A^^ is independent of the number x of 
linking branches (or the number m — x of free branches) . 

Matrices and vectors. One finds 



A. MFPT averaged over starting sites of level k 

The target is set at a*^^^. Let us introduce the param- 
eter z^i^, that numbers the repetition of a branch i^ : 

I'l = 1^2 = ^ — X 
^3 = ^4 = ^5 = ^ 

Due to the branch repetitions, it is necessary to adapt 
equation (23) by using the following notation : 



Mtot = ^^qMq ; Vtot =^^q^q ^tot = ^ J^q- 
q=l q=l q=l 



The analog of equation ( 23 ) is then : 

k 

E "^'"^ = ^0 ^tot H„ + J2n vUm\-1 Vtot , (33) 

paths l — l 



with Vi, 



. Let us define: 



/i = 3 



6m 



Ai = X, A2 = 2m + X 



1 



A^i = ( 1 Q 1 , -^2 = ( Q ^ ) , Ms 



1 



1 
2/3 1/3/ ' 



vi=i;),v.=(;),v,=(»),v.=i-i,v. 



Crossing time. One finds 



6m\ 



n—k 



(30) 



2^ 
(31) 



(32) 



where Ai and A2 are the two eigenvalues of Mtot- Thus: 

k 



paths 



+ ^2,1 



1^1 



A^-^(m -x)\^ J + A^-^(m + 3x) ^ 
1 



m — X . 
2V^ + ^)'^ 



1 



1 . + 3a; ^ ^fc / 1 



2 (/x-l)A2, 
m — a; /— 1\ m + 3x 



2X2 



1 



(m-i)A2 



-3a: /1\ \ /xxfc 



(34) 
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In order to relate the previous formula to (T) , we 
need to take into account the repetitions of branch z/g = 5, 
otherwise the nodes belonging to the levels < k — 1 would 
be counted several times. Let us call the number of 
sites belonging to level k. An examination of the labeling 

scheme shows that (T) is related to the 2nd coordinate 



{N, - i)lT)''^ = y: rc^) - - 1) E E -T^'' 



k-1 



paths 



1=0 paths 



which yields: 
• for X ^ 3: 



(35) 



(T) 



(fe) 



II" 



iVfc-1 



1 / m — X 



1 m + Sx \ A$+^ - xAS + a; - 1 



(m - 1)A2 



Ao - 1 



m — X m + 3x \ 3x — 3 — 2x {x/3) 



2A2 



(m - 1)A2 



X — 3 



(36) 



• for X = 3: in the latter expression the fraction 

(in the third term of the bracket) 
needs to be replaced by 1 — 2k. 



In expression ( 36 ) the volume of level k writes : 



Nk 2m 



A 



2. 



(37) 



B. MFPT averaged over the starting sites of 
connectivity 1. 

Let us keep the target at a'^^^ and calculate the MFPT 
averaged over the starting sites of connectivity equal to 
1. They belong to level n; let us call Nf^=i their number. 
One has: 



A^^=i =2{m-x) X^- 



(38) 



An examination of the labeling scheme leads to tlie_ fol- 
lowing formula, that relates the sum (A/'^=i — 1) (T) ^^-^ 
to the coordinates of J2paths '^^^^ ' which has been calcu- 
lated previously: 



l)(T>„=i= E (^'"^ 



paths 



E (mr("-i)| - (2a; - m) r'"-^) I ). 

(39) 



paths 



The MFPT averaged over all sites of connectivity 1 can 
then be calculated, using expression (34): 



(iV.=i - 1) (T),^, = 



3m + 2a; 



3m 



7x) 



X 



■x) 



A2(m 



!)■ 

(40) 



We show in appendix that this method also applies to 
the case of the (li, 'L')-flower networks introduced in [25] 
as examples of deterministic scale- free networks, that are 
either fractal or small-world. 



III. A DIFFERENT METHOD FOR 
HIERARCHICAL GRAPHS 

A. Recursivity and Modularity 

We now consider a different class of graphs, that is hier- 
archical, non-decimable, self-similar networks, which are 
built deterministically and recursively in a manner remi- 
niscent of exact fractal lattices. More precisely, the graph 
of generation g is obtained by properly Unking together 
a certain number of copies of generation g — 1. Differ- 
ently from networks previously analyzed, where different 
replicas meet at a single node, here exact renormalization 
procedures are not applicable. Yet, we can exploit modu- 
larity to detect analogous subgraphs whose nodes satisfy 
intrinsic, mutual relations, and self-similarity, which al- 
lows to establish recursion relations. 

Now, in order to fix ideas we focus on a particular 
example of hierarchical network, introduced in [26] and 
further investigated in [27- 30 , (see Fig. [7|; other exam- 
ples can be found in ^22^ .31 -.34 . 

By denoting as Qg the graph of generation g, we have 
that Qo is given by a single node, also called "root" , while 
Qi is a chain of length three obtained from Qq by adding 
two more nodes and connecting each of them to the root; 
the two nodes added are called "rims" (of level 1). Sim- 
ilarly, at the second iteration, one introduces two copies 
of ^1, whose rims are directly connected to the root: now 
the root is connected to the original two rims of level 1 
and to four rims of level 2. 

Proceeding analogously, at the ^-th iteration one in- 
troduces two replica of the existing graph, i.e. of Gg-i, 
and connects the root with all the new 2 x 2^~^ rims, 
referred to as rims of level g. Hence, the root turns out 
to be a hub connected with 2'^ rims of level n, where 
n G [1,^], in such a way that its coordination number is 
Zg = 2(2^ — 1), on the other hand, rims of level n have a 
coordination number equal to n. 

Given a rim of level n, here referred to as r^, one can 
see that it is not only connected to the root but also 
to other "minor hubs" hk,n^ namely nodes that work as 
main hub for any subgraph Qi^^ k = l,...,n, containing 
both Tn and hk^n] more precisely, we refer to hk^n as the 
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hub of height /c, with respect to a rim of level n (see 
Fig.[7|). The root will be also referred to as the main hub 
and denoted as = hg^g. Also, given a node i which is 
a rim of level n, we say that the set of rims of the same 
level and sharing with i the same hub of height k are rims 
shifted by k with respect to i; this set is denoted as {vk^n} 
and its cardinality is |{r/e,n}| = 2^ (see Fig.[7|). The total 
number of nodes making up Qg is Ng = 3^, while the 
total number of rims is ^^^^^ 2^ = 2(2^ -1) = Zg. 

Furthermore, we mention that the degree distribution 
for hubs is given by the power law P{k) ^ k~^ ^ with ex- 
ponent 7 = log 3/ log 2 ^ 1.59, while the remaining nodes 
follow an exponential degree distribution P{k) ~ (2/3)^. 
For further details about the topological properties of Qg 
we refer to [26-30 . 

We also notice that each subgraph making up the 
whole graph can be looked at as a module; connections 
between a module and the remaining graph are few (with 
respect to the size of the subgraph itself) and concern 
only the pertaining rims. Indeed, it is possible to de- 
termine a hierarchy of nodes, based on their degree of 
clustering, consistently with [22]: Although for Qg it is 
not possible to establish a one-to-one correspondence be- 
tween the clustering coefficient of a vertex and its degree, 
one can see that the clustering coefficient [38 decreases 
with the degree. 

The MFPT's on Qg have already been analyzed for 
special target locations [ 2M3Q] and, before proceeding, 
it is worth recalling some results which may be useful in 
the following. In particular, for a simple RW on Qg^ the 
mean time to ffist reach the main hub H starting from 
an arbitrary rim of level g is [28 

T,(F,rg) = ^Qy-3, (41) 

while the mean time to ffist reach any of the 2^ rims of 
level g starting from the main hub is 

T,{{r,},H)=^-(^^' -I, (42) 

where the mean is taken over all possible paths; no- 
tice that the asymptotic behavior of Tg{H^rg) and 
Tg{{rg}^ H) is the same, namely ~ (3/2)^, even if the 
number of targets is 1 and 2^, respectively. 

B. Labeling code and time to main hub 

We now introduce a method which allows to calculate 
straightforwardly the mean time to ffist reach the main 
hub starting from an arbitrary node. First of all, we need 
a proper labeling for nodes which exploits the topological 
features of the structure. Basically, we associate to an 
arbitrary node i belonging to the graph Qg a code, e.g. 

= {lrt..rrt)^ made up of g letters properly chosen in the 
alphabet {t,r, as we are going to explain. The whole 
graph can be looked at as the combination of three graphs 



of the previous generation: Qg-i (corresponding to t) and 
two copies of Qg-i arranged on the right (r) and on the 
left (/), respectively. Now, according to which of these 
main subgraphs i belongs to, we have that is equal to 
either t, r or Once detected the main subgraph, one 
proceeds analogously distinguishing the three subgraphs 
of second order, i.e. Qg-2^ and evaluating which contains 
the node hence determining ^f. Finally, at the ^-th 
iteration, one is left with the subgraph ^i, in such a way 
that its subgraphs are simply three nodes, one of them 
corresponds to i. For instance, referring to Fig. [7| we 
have: 

68 = {ittt), 

= m). 

^2 = {Irrt). 

As anticipated, we focus on arrangements where source 
and target belong to different modules; as we will see, 
this typically requires the passage through in such 
a way that we ffist need to calculate the time T{H^i). 
For this aim, our addressing, while able to determine 
univocally a node, is somehow redundant, since, due to 
the intrinsic symmetry, the distinction between left and 
right subgraphs is unnecessary. For this reason, one can 
denote any of the two subgraphs in the bottom as b in 
such a way that, for a graph of generation ^, one could 
write 6 = {t'^b'\..t'^-^b'^), with Yl^^i k = 9] also notice 
that, without loss of generality, we can always assume 
(^l = t and = 6, with ii and ik possibly zero, while 
ii > for I e [2, /c — 1]. Reading this string from right 
to left, we can write a general expression for the MFPT 
from i to H. In fact, assuming that z is a rim {^f = b) 
of a certain inner subgraph Qg^ , in order to reach we 
need to pass through the main hub of Qg^ itself, where gi 
is simply ik- Now, the main hub of Qg.^ is either H (when 
ik = g or when /c = 2) or the main hub of a certain inner 
subgraph Qg^^ where g2 turns out to be ik + ik-i- One 
can proceed analogously, bouncing from hub to rim and 
from rim to hub of larger and larger subgraphs, in such a 
way that the following general expression for the MFPT 
from i to H holds 

T{H,i) = T{H,ii)= (43) 

k 

= J2[Tj,{H,rj,)+Tj,^,{{rj,^,},H)], 

1 = 1 

I 

31 =9 - ^k'^ 

l' = l 

where we used H to indicate the main hub of the 
(sub)graph considered (denoted by the index ji) and 
To = 0. Recalling the examples above and using Eqs. [4l] 
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FIG. 7: (Color on line) The graph Qg of generation g — 4:. Darker nodes are hubs and brighter nodes are rims (of different 
subgraphs/levels). The labeling is complete for the subgraph in the top, and proceed analogously for the other subgraphs for 
which only a few labels have been inserted as example. Here 1 represents the main hub, i.e. H = /i4,4, nodes {2,3} are rims 
of level 1, nodes {5, 6, 8, 9} are rims of level 2 and so on, as specified. Also, focusing on node 14, we notice that 10 and 13 are 
/i2,3 and /ii,3, respectively, while node 15 is ri,3, {17, 18} = {^2,3} and {23,24,26,27} = {^3,3}- 



and [42} we write 

T{H,U) = T{tHt) = m{n},H)+T2{H,r2), 

T{H, ^28) = T{bt^) = Ts{{rs},H)+Ti{H,ri), 

T{H,^52) = T{bH) = T^{{n},H)+n{H,r^). 



Of course, summing up such times over the whole set of 
nodes, one recovers the global mean first passage time 
ra = Er^HnH,i)/{N-l) EH EI]. 

Now, in order to complete the calculation for the 
MFPT from an arbitrary source z to a target j, we need 
T(j, which can be calculated exploiting the centrality 
of H and the self-similarity of the graph, by implement- 
ing a set of recursive equations. In order to preserve the 
generality of the method we focus on a particular class of 
targets, easily identifiable in generic hierarchical graphs, 
namely on rims of an arbitrary level n. 



C. MFPT's from hubs 

Beyond those discussed before, in order to calculate the 
MFPT from H to a rim r^, we need further quantities. 
First, let us introduce the following: Tg{rn^H), which 
represents the mean time to go from the main hub to a 
rim of level n, Tg{rn,hk,n)^ which represents the mean 
time to go to a rim of level n from a hub of height k 
with respect to r^, and Tg{rn^ ^k,n)i which represents the 
mean time to reach a rim of level n from a rim of the 
same level, but "shifted" by k (see Fig. [t]) [39]. 

Then, we can write the set of equations: 

1 1 """"^ 

Tg{rn^H) = - + -^2^[l + T,(r,,n+i,0] (44) 

^ n — 1 

^3 1=1 
1 ^ 

^9 l=n+i 
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where the first term in the r.h.s. accounts for a direct 
jump from the root to the target, the second one accounts 
for shifted rims of level n itself and the remaining terms 
account for rims of all levels other than n; 



k-l 



Tg{rn,hk^n) = - + -^2^[l + T^(rn,n+i,n)] (45) 

^ k-l 

+ -^2'[l^Ti{H,ri)^Tg{rn,hk,n)]. 
1=1 

similarly to the previous case; 
1 

Tg{rn,rk,n) = -^[l^Tg{rn,hi^n)] (46) 
^ i=k 

k-l 

+ -Y,[l^Tg{rn,rk,n)^Ti{{ri},H)], 

^ 1=0 

where the first term in the r.h.s. accounts for the fact 
that to reach you need to pass through a common 
minor hub, or, possibly H itself, while the second term 
accounts for bounces from the starting point to close (non 
common) minor hubs. 

This system of recurrent equations can be solved 
by first focusing on Eqs. [45p6| and building up the 
differences between terms for k -\- 1 and k so to get 
rid of the sums. The solutions found for Tg{rn-, hk^n) 
and Tg{rn-,rk^n) are then plugged into Eq. 44 and, ex- 
ploiting Eqs. 4T][42l as well as proper initial conditions 
(e.g. 2T(rn,/ii,n) = 2 + T(rn,ri,n) and Tg{rn,n,rn) = 
Tg{H^rn) + Tg{rn-,H))^ one obtains closed form expres- 
sions which read as: 



1 

n-1 



-V^(n) 



•2^-^(3^-2^) 



(47) 



-6{k,n) (48) 



1 



n 



n 



(l){k,n) 



, (49) 



2^ n-k ' 2 
where 

Cn = 8(3^ - 2^)[2n - 1 - 2^-^n0(n - l,n) + nV^(n)]-\ 
and 

k-l 



(t)(k,n) = 



i=0 



n — 1 



n-l . 
1=0 ^ 



(t)(l,n] 



All these formula have been successfully checked ver- 
sus numerical estimates obtained by means of the pseudo 
Laplacian ^35j . 



It is convenient to report the asymptotic (A^ ^ oo, i.e. 
g large) expressions of previous quantities, which turn 
out to be the same for all of them, namely: 



3^ 



(50) 



so that it is easy to see that the level n plays algebraically: 
although the distance hub-rim and rim-rim is equal to 1 
and 2, respectively, whatever n, /c, rims added at larger 
generations are "easier" to be reached. We also notice 
that the height of the minor hub considered or the shift 
among rims, just provide minor order corrections. In par- 
ticular, once generation and level are fixed, the MFPT 
decreases with /c, the reason is that, although the distance 
between starting point and target remains the same, in- 
dependently of /c, a large k implies the passage through 
more connected hubs which are easier to be reached. 

Analogous recursive equations can be implemented for 
the case of multiple targets (see for examples ^28| |29 ] ), 
while here we just focus on the case of single target. 



D. Examples 

The results explained in the previous section, together 
with those summarized in the Sec.VA, allow to get an 
exact expression for the MFPT between two nodes i and 
j, such that the path has to include a hub h. In this way 
one first detects the hub h and then calculate T{i^j) as 
a sum of the partial MFPT from i to h and from h to j. 
In order to clarify the procedure, we now present some 
examples where several kinds of situations are considered. 



Source and Targets are both rims 

Being i a rim of level n and j a rim of level m with n ^ m 
[W, it is easy to see that, in order to go from i to j (or 
vice versa), one has to pass through so that 

T(j,z) = T{rm.rn)=Tg{H,rn)^Tg{rm.H) 



Tn{H,rn)^Tg{r^,H) 



1 

m—1 



m 



3 V2 

+ 3 



22-m/3^ - 2^) - — . 

m 



(51) 



Hence, the leading term is typically Tgirm-, H). This also 
implies that the MFPT for the same nodes, but opposite 
direction, i.e. T{rn,rm)^ differs from T{rm^Tn) and their 
ratio goes like n/m. 

For example, let us refer to Fig. [7| and let us fix i = 14 
and j = 41 (or, of course, equivalent nodes). Then, we 
find 

T4(41,14) = Ts{H,rs)^T4{r4,H) = 6+1109/12 ^ 98.42, 
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where we used Eqs. |4T]and|47l Analogously, 

T4(14,41) =T^{H,u)^T^{r3,H) = 21/2+809/6 ^ 145.33. 

As expected from Eq. |5lj T(14,41) > T(41, 14) due to 
the fact that from the main hub it is easier to reach a 
rim which belongs to higher levels. 

Source and Targets are a hub and a rim 
Being i a hub of height k with respect to a set of rims 
of level n and j a rim of level m with n ^ again, in 
order to go from i to j (or vice versa), one has to pass 
through so that 



T{j,i) 



Tg{{rn},hk,n) 

Tk{{rk},H)^ 



+ Tg{H,r„)+Tg{rm,H) 
Tg{H,r„)+Tg{rra,H) (52) 













-1 + 


1 


1 


m 












-2^) 



'tp{m) 
m ' 



3-4 



(53) 



Hence, the leading term is typically Tg{rm^H). This is 
rather intuitive as Tk{{rk}^H) accounts for trapping on 
a set of 2^ nodes, Tg{H,rn) for trapping on the main 
hub, while Tg{rm-,H) for trapping on a single node with 
relatively small (equal to m) coordination number. 

As an example, let us consider the graph ^4 in Fig. [7| 
and calculate 

7 21 809 
T4(14,28) = 2 + y + "6" ^ 

For comparison, let us also consider for a "rescaled 
case", that is 

T3(5,10) = 2 + 6 + 54 = 62. 



Source and Targets are rims of the same subgraph 
Let us consider the subgraph Qg' of Qg^ where g' < g 
and let us fix i and j as rims of level m < n < g' oi the 
subgraph, whose main hub has to be crossed in order to 
go from i to j; for simplicity, let us assume that j also 
belongs to the set {r^} in Qg. Therefore we can write 



T{j,i) 



Tg,{H,r„ 

8 /3^ 

3 V2- 
3^ 



Tg{rn^ hg' 
Tg{rn^ hg' 



(54) 
(55) 



-(t){k,n) 



For example, still referring to the labeling of Fig. [7| let 
us fix i = 32 and j = 41. Then, for ^ = 4 we find 

T(41,32) = T2(i^,r2)+T4(r4,/i3,4) = 3+887/12 ^ 76.92, 

where we used Eqs. [4l] and |48) Analogously, for ^ = 3 
we find 

T(14, 11) = Ti(i7, n) + T3(r3, /i2,3) = 1 + 211/6 « 38.17. 



IV. CONCLUSIONS 

In this work we introduced general methods to calcu- 
late exactly first-passage quantities on self-similar net- 
works, which are defined recursively. In particular, we 
focused on the mean first-passage time from a source S 
to one or to several target sites {T^}, and on the split- 
ting probability, namely the probability to reach, starting 
from a given target before all other targets. Indeed, 
these quantities allow a sound description of a wide range 
of dynamical processes such as diffusion limited reactions 
or search processes embedded in complex media [36 ,^. 

In general, the methods introduced strongly rely on 
the recursivity of the underlying structure, namely on 
the fact that the whole graph can be built according to 
a recursive procedure: at the n-th iteration the graph 
Qn is obtained by properly combining a finite number of 
graphs Qn-ii each corresponding to the graph itself at 
the previous iteration. 

Recursive networks where replicas meet at a single 
node are amenable to exact analysis by renormalization 
techniques and we show that the above-mentioned first- 
passage quantities can be recast as solutions of simple 
matricial equations. We also consider examples where 
replicas are connected by links and exact decimation is 
no longer accomplishable; then, one can impose a num- 
ber of coupled equations, each corresponding to a set of 
equivalent nodes. 

Hence, a broad range of topologies can be addressed 
via a unifying approach: we considered explicitly recur- 
sive networks as diverse as finitely ramified fractals (Sier- 
pinski gasket, T-fractal), scale-free (trans)fractals {{u.,v)- 
fiowers), non-fractals, mixtures between fractals and non- 
fractals (Song-Havlin-Makse networks), non-decimable 
hierarchical graphs (Barabasi-Ravasz-Vicsek network). 

In any case, calculations performed are exact; results 
previously obtained for special cases of target arrange- 
ments are recovered and extended to account for more 
general configurations. 



Appendix A: The T-graph 

In this section we apply the method that has been pre- 
sented in section [l] to the T-graph (see Fig. [T]). We first 
give the matrices M.i^^ the vectors Vi^ and the exit time 
T/c. It is then possible to calculate explicitly splitting 
probabilities and MFPTs for any starting point on the 
network. For example, we calculate explicitly two differ- 
ent types of averages: the MFPT averaged over starting 
points of a given level and the MFPT averaged over all 
starting sites with a given connectivity. In the particu- 
lar case where the average is performed over all starting 
positions, we recover the result of Ref ^10] . 
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^(fc-i) 



ik = 2 



ik = S 



FIG. 8: Renormalization scheme and indexes for the T- graph 



General results 



Quantities V and T are here 2-dimensional vectors, 



and will be denoted 



(k) 



Figure 



(but not unique) labeling scheme {p = 
Matrices and vectors. One finds 



shows the chosen 
3). 



Ml = 



1 
1/2 1/2 



Mo 



_ (1/2 1/2 



M. 



_ A/2 1/2' 
1/2 1/2 



(Al) 



As a result of the indexation scheme that has been 
chosen, S^^^ = Y^paths ^^^^ 5 ^^^^^ words the MFPT 
summed over all sources of level k is equal to the 2nd 
coordinate of Eq. (A5). Therefore, using the fact that 
there are 3^ sources, one gets : 



(T) 



7 



6x3^ 2 X 6^+1 



(A6) 



For the particular case /c = n, the latter expression gives 
the exact result derived in [lO] . 

MFPT averaged over starting points of given con- 
nectivity. The nodes of the T-graph hav e connectivity 
Ai: = 1 or 3. It is possible to use (A5) to determine 
the MFPT averaged over all the N,^=^ starting sites of 
connectivity n = ?> (which of course gives access to the 
similar quantity for the starting sites of connectivity 1). 
Indeed, an examination of the labeling scheme applied 
to th e entire network shows that the first coordinate of 



(|A5|) is equal to 2 N,^^^ (T)^^^. Then, using the fact that 

3 = — ; 



one obtains : 



1 



1 



2 X 3^+1 1 



Crossing time. One finds 



V3 



(A2) 



(A3) 



2. Examples of averaged MFPTs 

Let us consider the case where the target is located 
at the extreme left vertex of the T-graph : T = a^^^. 
Note that an average MFPT calculated for that target 
for a network of generation n is equal to an average for 
the next (n + 1) generation with a target at the center. 
Using that setting, the equivalent of equation ( 15 ) is: 



rC') = ^ T, {M,, ■ ■ ■ M,,^,) V,, with Vi, 



1=0 



(A4) 

MFPT averaged over the sources of level k. Follow- 
ing the method developed in section [iDj we make use of 
Eq. (23) and obtain: 



paths 



ik) 



[4-2-1 [i 

4 X 3^ - 3 X 2-^^+^) 



(A5) 



Appendix B: The (u,v)-f[owers 

In this section, we consider the case of the (li, v)-^owei 
networks introduced in [25] as examples of deterministic 
scale-free networks, that are either fractal or small-world 
depending on the values of the two parameters u and v. 
The algorithm to build the {u^ i')-flower networks is as 
follows: one starts (generation n = 0) with two sites con- 
nected by a link; generation n + 1 is obtained recursively 
by replacing each link by two parallel paths of respec- 
tively u and V links. In order to simplify some notations, 
we define w = u-\-v^ and take (without loss of generality) 
u <v. Examples of (1,3) and (2, 2) -flowers are shown in 
Figure |9] Let us remark that the connectivity of a site 
is determined by its level: = 2^~^~^^ when k > 1 and 
^0 = 2". 

We focus on the examples {u = l,v = 3), {u = 2^v = 
2) and {u = 2^v = 3). As will be shown in Section B 1 



the (l,3)-flower is a transfractal (small-world) network, 
whereas the (1,3) and (2, 3) -flowers are fractal networks. 
Note that the GMFPT for a target on level was previ- 
ously obtained by Zhang et al. [12] in the cases (1,3) and 
(2, 2) flowers. We show here that the MFPT can be ob- 
tained for any {u, 'y)-flower network for any starting site, 
and that averages over different classes of starting sites 
can be obtained. 
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(a) 



• • y 4 ^ 



D 




C 








n=2 






' ' A 




B 











n=3 








r 



^(k)i l5(kyk)i 1 i 1 i ^k) 

«k = ik = 1 «k = 2 ik = 3 

(i) M - 1 , w - 3. 

^(k-i) ^(k-i) ^(k) y\ /\pW 



«k = ik = 1 

(ii) M = 2 , - 2. 



«(k-i) ^(k-i) 




ik = ik = 1 ik = 2 

,5(k) 

ik = 3 ik = 4 

(iii) M - 2 , - 3. 



FIG. 9: Renormalization scheme for two (u,v)-B.ower net- 
works : (a) u = l,v = 3 ; (b) u = 2,v = 2. Figure re- 
produced from [19]. 



FIG. 10: Color online, (i^, t')- flower networks : indexes. 



2. (l,3)-flower 



1. General properties of the (u^v) flowers 

The number of links of level k is and the number 
of sites of level k is given by |25] : 



w-2 
w — 1 



w 



(Bl) 



The value of the parameter u splits the scaling of the 
diameter R with respect to n into two distinct cases: 



R' 



{v — l)n for 11 = 1, 
u'^ for u> 2. 



(B2) 



Therefore, when = 1, the network is a small- world. 
When > 2, it is a fractal network, and its fractal di- 



mension is obtained by combining equations (Bl) and 
(|B2|): 



hi{u + v) 
Ini^ 



(B3) 



Crossing time. Following the method presented in 
Section [TCl we find a general formula for r/^: 



n — k 



u-\- V 



2 u—l v — l 

This result holds also for = 1. 



(B4) 



We first study the case = 1, = 3, as an example 
of small world network. Figure 10 (i) shows the labeling 
scheme that has been chosen. We find that Tk = 3^~^ 
and: 



Mo 



M2 



1 
1 

2/3 1/3 
1/3 2/3 



Vo 



Ml 
Ms 

V2 = 



1 

2/3 1/3 

1/3 2/3 
1 

, V3 = 



(B5) 



(B6) 



It is then possible to calculate exactly the splitting prob- 
abilities and MFPTs for any starting site on the network. 

MFPT averaged over starting sites of level k. We as- 
sume that the target is located at a^^^ in leve l 0, and we 
apply the method d evelo ped in Sections |ID| (derivation 
of ^paths T^^^) and II A (substraction of branches repe- 





titions). We first use formula m\ with T^^^ 
derive 

' -1 



paths 



1 



1 



1 



2x3^ 



3n 



to 



(B7) 

We then focus on the 2nd coordinate of the latter expres- 
sion, but we need to take care of the fact that branches 
i/c = and i/c = 3 generate two contributions of the same 
source a way to avoid this repetition is to substract 
the same quantity for all the levels lower than /c, just like 
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in Eq. (ISSl) (with x = 1). We then get: 



the chosen labehng scheme. One finds 



(Nk - 1) (T) 



2x3n-l4/c^3n-AC4/c_2x3^ 



k Ak 



(B8) 



with Nk = I (4^ + 2) . For the particular case k = the 
latter expression is in agreement with [12 . 



3. (2,2)-flower 

In this example we find =4"^"^ and: 

= (l/2 1/2) ' 



1/2 1/2 
1 



(B9) 
(BIO) 



It is then possible to calculate exactly the splitting prob- 
abilities and MFPTs for any starting site on the network. 
MFPT averaged over starting sites of level k. Again, 

the target is located on a*^^^ and T*^^^ = ^4^^ • "^^^ 

repeted branches must be taken into account: the 
method has been given in Section II A Here uq = ui = 2. 
Using formula (33) we get: 



paths 



-)k-l 



-1 



6 



4^+1 



1 



(Bll) 

Then we focus on the 2nd coordinate of the latter expres- 
sion, and apply equation (35): 



(Nk - 1) (T) 



(fe) 
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[2 X 4'^+! + 3fc + 10] (B12) 



with Nk = ^ (4^^ + 2) . For the particular case k = n, the 
latter expression is in agreement with [12 . 



•^0- 1^2/3 l/sj' -^1- 1^1/3 2/3; ' 
=- 2/3 



1 / ' 

1 



(B13) 



■^'=[1/2 1/2 j' ^'={%^ /' 



(B14) 



It is then possible to calculate exactly the splitting prob- 
abilities and MFPTs for any starting site on the network. 

MFPT averaged over starting sites of level k. The 
target is located on a*^^^. We again apply formula (23) 

with r^^) = 1 and find : 



paths L \ / \ / \ 

(B15) 

Substracting the same quantity for all the levels lower 
than k (see Eq. (|35|) finally yields: 



(Nk - 1) (T) 
with Nk 



(k) 



- 5'=-! - 7 X 5'=-! e-'^ + - 



(B16) 



3x5^+5 
4 
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